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ABSTRACT 

Starburst99 is a comprehensive set of model predictions for 
spectrophotometry and related properties of galaxies with active star 
formation. The models are an improved and extended version of the data 
set previously published by Leitherer & Heckman (1995). We have upgraded 
our code by implementing the latest set of stellar evolution models of the 
Geneva group and the model atmosphere grid compiled by Lejeune et 
al. (1997). Several predictions which were not included in the previous 
publication are shown here for the first time. The models are presented in 
a homogeneous way for five metallicities between Z = 0.040 and 0.001 and 
three choices of the initial mass function. The age coverage is 10 6 to 10 9 yr. 
We also show the spectral energy distributions which are used to compute 
colors and other quantities. The full data set is available for retrieval at 



http : //www . stsci . edu/science/starburst99/| . This website allows users to 



run specific models with non-standard parameters as well. We also make the 
source code available to the community. 

Subject headings: galaxies: evolution -- galaxies: fundamental parameters - 
galaxies: starburst — galaxies: stellar content 



-4- 



1. Introduction 

Synthesis models for comparison with observed galaxy properties, and in particular 
with spectral energy distributions, have become increasingly popular in recent years. 
This is partly due to improvements in stellar libraries and models. An overview over 
the latest achievements can be found in the conference volume edited by Leitherer, 
Fritze-v. Alvensleben, & Huchra (1996b). At least equally important is the possibility 
of electronic dissemination of the data, either by CD-ROM or via the internet. Most 
groups who provide such models to the community make them available electronically. See 
Leitherer et al. (1996a) for a compilation of a suite of models. 

This paper concerns a particular set of synthesis models: models which are optimized 
to reproduce properties of galaxies with active star formation. In the absence of an 
active galactic nucleus, most radiative properties of such galaxies are determined by their 
massive-star content. The most extreme examples are referred to as starbursts (Weedman 
et al. 1981; Weedman 1987; Moorwood 1996) but the models in this paper are often 
applicable to less extreme star-forming regions and galaxies, like 30 Doradus and M33, 
as well. Leitherer & Heckman (1995; hereafter LH95) published a homogeneous grid of 
synthetic starburst models. These models, together with those of Bressan et al. and 
Bruzual & Chariot (both published on the CD-ROM of Leitherer et al. 1996a), Fioc & 
Rocca-Volmerange (1997), and Cerviho & Mas-Hesse (1994), are widely used to aide the 
interpretation of galaxy observations. 

Although the LH95 models are still recently up-to-date and no major error was found, 
we decided to generate and release a new generation of models. There are two reasons for 
this. First, stellar modeling is continuously advancing, and some known shortcomings of the 
stellar atmosphere and evolution models used in LH95 have been improved. Specifically, 
we have included the new model atmosphere grid of Lejeune, Buser, & Cuisinier (1997) 
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and the latest Geneva evolutionary tracks in our code. Other modifications include the 
possibility to perform isochrone synthesis, a technique which was pioneered by Chariot & 
Bruzual (1991). 

The second motivation for an update to LH95 concerns the distribution method. 
While the LH95 data set is unofficially available via our website, the principal publication 
is on paper. Many users found this inconvenient and encouraged us to publish fully 
electronically. Therefore we decided not to include any tables with model predictions and 
only a subset of the figures in the hardcopy version of this paper. The full data set with 
all the figures is available in the electronic version of the journal and via our webpage, 
which is an integral part of this publication. The reader should visit the website at 
http : //www . stsci . edu/science/starburst99/| to view the figures and to download the 



tables. Moreover, this website offers the opportunity to run tailored models remotely and 
to access our source code. The reader should use this paper as a guide when visiting our 
webpage. 

The organization of the paper is as follows: In § 2 we outline our computation 
technique, the model assumptions, and the improvements over LH95. This section also 
explains how to reach our website and how to navigate. The spectral energy distributions 
are presented in § 3. Spectral line profiles in the ultraviolet (UV) are discussed in § 4. In 
§ 5 we provide numbers for the stellar inventory. Luminosities and colors are in § 6 and 
7, respectively. Properties related to the far-UV are presented in § 8. Some other useful 
diagnostic lines are in § 9. § 10 covers the mass and energy return of a stellar population. 
We conclude with 5 11. 
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2. Model assumptions and computational technique 

The present model set is an extension of our previous work which was published in 
several papers. Most of the earlier tables are in LH95. The energy distributions used in 
LH95 were made available electronically in Leitherer et al. (1996a) and on our website. 
Ultraviolet spectra between 1200 A and 1850 A at 0.75 A resolution were originally published 
in Leitherer, Robert, & Heckman (1995) and distributed electronically in Leitherer et al. 
(1996a). A set of spectra around O VI A1035 and Ly/3 was discussed by Gonzalez Delgado, 
Leitherer, & Heckman (1997). All mentioned models have been re-computed and included 
in the present database. The covered parameter space is similar to the one in LH95. Since 
arguments are given in LH95 for the choice of these particular parameters, we will be brief 
and only address those issues in more detail which are different from LH95. 

We consider two limiting cases for the star-formation law: an instantaneous burst 
of star formation and star formation proceeding continuously at a constant rate. The 
instantaneous models (also referred to as a single stellar population) are normalized to a 
total mass of 10 6 M . The star-formation rate of the continuous model is 1 M Q yr _1 . These 
normalizations were chosen to produce properties which are typical for actual star-forming 
regions in galaxies. 

As we did in LH95, we offer three choices for the stellar initial mass function (IMF). 
The reference model is a power law with exponent a = 2.35 between low-mass and high-mass 
cut-off masses of M\ ow = 1 M and M up = 100 M , respectively. This approximates the 
classical Salpeter (1955) IMF. Most observations of star-forming and starburst regions are 
consistent with a Salpeter IMF - - although the uncertainties can be large (Leitherer 1998; 
Scalo 1998). For comparison, we also show the results for an IMF with a = 3.3 between 1 
and 100 M , which has a higher proportion of low-mass stars. Finally, we computed models 
with a truncated Salpeter IMF: a = 2.35, M\ ow = 1 M , M up = 30 M . The models can 
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easily be re-normalized to other M\ ow values. The lower mass cut-off is often only a scaling 
factor since low-mass stars have a negligible contribution to the properties of the stellar 
population, except for the mass locked into stars. Absolute quantities (such as, e.g., the 
number of ionizing photons) can be rescaled to other M\ ow values by multiplying our results 
by factors of 0.39, 1.00, 1.80, 3.24 for M low = 0.1, 1.0, 3.5, 10 M . Relative quantities (like 
colors) should, of course, not be scaled. Note that 10 M stars are no longer negligible for 
many model predictions, and the models should be recalculated for this cut-off mass, rather 
than scaling the existing data. Common sense is necessary when applying this scaling factor 
to instantaneous models in particular. The lifetime of a 2 M Q star is only somewhat above 
1 Gyr. Therefore burst models at 1 Gyr are actually only represented by a population of 
1-2 M stars. Clearly even a small change of the lower mass cut-off will drastically alter 
the model predictions. 

We have implemented the new set of stellar evolution models of the Geneva group - 
as opposed to the Maeder (1990) models in LH95. For masses of 12 -25 M Q (depending on 
metallicity) and above we are using the tracks with enhanced mass loss of Meynet et al. 
(1994). The tracks of Schaller et al. (1992), Schaerer et al. (1993a, 1993b), and Charbonnel 
et al. (1993) with standard mass loss are used between 12 and 0.8 M Q . We have favored the 
enhanced mass-loss models over the standard ones because they are a better representation 
of most Wolf-Rayet (WR) star properties, except for the mass-loss rate (M) itself. The 
high mass-loss rates of WR stars are in conflict with observations (Leitherer, Chapman, & 
Koribalski 1997). This may suggest the lack of one or another ingredient in the evolution 
models, like rotation induced mixing processes (Maeder 1995). The evolutionary tracks 
do not take into account binary evolution. Inclusion of Roche-lobe overflow in binary 



2 A typographical error was introduced in eq. (3) of LH95: the numerator should be 0.80 
instead of -0.35. 



systems can modify some of the predictions, in particular during the WR phase. Models 
taking these effects into account were computed by Cerviho, Mas-Hesse, & Kunth (1997), 
Vanbeveren, Van Bever, & de Donder (1997), Schaerer & Vacca (1998), and Dionne (1999). 

The Geneva models include the early asymptotic-giant-branch (AGB) evolution until 
the first thermal pulse for masses > 1.7 M . Stellar evolution models in the AGB phase are 
rather dependent on input assumptions, and therefore quite uncertain. A recent comparison 
between different sets of AGB models was made by Girardi & Bertelli (1998). 

Five metallicities are available: Z = 0.040, 0.020 (= Z ), 0.008, 0.004, and 0.001. 
These are the metallicities of the evolutionary models. All element and isotope ratios are 
independent of metallicity. Note that we do not treat chemical evolution self-consistently: 
each stellar generation has the same metallicity during the evolution of the population. The 
error introduced by this simplification is negligible as long as star formation proceeds for 
less than about 1 Gyr, otherwise self-consistent models are needed, such as those by Moller, 
Fritze-v. Alvensleben, & Fricke (1997). 

The models cover the age range 10 6 to 10 9 yr. In principle we could have evolved the 
populations up to a Hubble time but we decided not to do so. The evolutionary models are 
optimized for massive stars and become less reliable for older ages. Many properties can 
still be safely predicted beyond 1 Gyr whereas others (such as, e.g., near-infrared colors) 
become increasingly uncertain. We will comment on uncertainties when the individual 
figures are discussed. 

Atmosphere models in LH95 were from Kurucz (1992) for most stars, supplemented by 
the Schmutz, Leitherer, & Gruenwald (1992) models for stars with strong winds, and black 
bodies for the coolest stars. Lejeune et al. (1997) have produced a new, homogeneous grid 
of atmospheres covering the entire Hertzsprung-Russell diagram (HRD) populated by the 
evolutionary models, including the coolest stars. This new grid has been implemented in 
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our synthesis code. We used the corrected grid whose fluxes produce colors in agreement 
with empirical color-temperature relations. The original atmosphere models do not include 
metallicities of 0.008 and 0.004. We interpolated between adjacent metallicities to obtain 
grids whose metallicities match those of the evolutionary models. No interpolation was 
required for 0.040, 0.020, and 0.001. As in LH95, the Schmutz et al. (1992) extended 
model atmospheres are used for stars with strong mass loss. The prescription for the switch 
between extended and plane-parallel atmospheres is the same as in LH95. The extended 
atmospheres are assigned to all stars with surface temperature >25,000 K and hydrogen 
content on the surface <0.4. We decided not to use extended model atmospheres for stars 
closer to the main-sequence since the benefit of improved wind treatment would be more 
than compensated by the limited metallicity range available for such models. 

The nebular continuum has been added to all spectrophotometric quantities, except 
to the energy distributions in Figures 7a-e, 8a-e, 9a-e, lOa-e, lla-e, and 12a-e. We have 
assumed that all the stellar far-UV photons below 912 A are converted into free- free and 
free-bound photons at longer wavelengths (case B). However, the stellar far-UV flux has not 
been removed in the spectral energy distributions, as this would probably not be desirable 
for the user. To alert the user, the hydrogen and helium emission coefficients of Ferland 
(1980) do not extend longward of 4.5 /xm. Therefore the nebular continuum at these 
wavelengths is undefined in our code. In practice, this is hardly relevant since dust emission 
often dominates at A > 3 /im (Carico et al. 1988) and our models become inapplicable 
anyway. 

All models shown here have been computed with the isochrone synthesis method, as 
opposed to the classical evolutionary synthesis method in LH95. Isochrone synthesis was 
introduced by Chariot & Bruzual (1991) as a technique to overcome the discrete appearance 
of the model predictions at late evolutionary stages when the resolution in mass of the 
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evolutionary models becomes inadequate. Instead of binning stars in mass and assigning 
them to a specific track, continuous isochrones are calculated by interpolating between 
the tracks in the HRD on a variable mass grid using a subroutine originally provided by 
G. Meynet. This scheme generally produces smooth output products, the only exception 
being quantities related to supernovae (see below). 

The time resolution of the model series is 0.1 Myr. Spectral energy distributions, 
however, are given only at time steps during which significant changes occur. 

The model predictions are discussed in the following sections. A summary of all the 
figures is in Table 1. Each figure has five panels (a-e) for each of the five metallicities. 
An exception are Figures 13 through 36 which show UV line spectra. These spectra are 
for Z only, as the library stars are currently available only for this metallicity. The 
models will be extended as soon as non-solar library stars become available (Robert 
1999). The data files which were used for the generation of the plots can be viewed and 
downloaded (see instructions on the website). The files with spectra (Figures 1 through 
36) contain data at more time steps than plotted. For clarity reasons we included only 
a subset of all time steps in the plots. Only the figures for solar metallicity models are 
reproduced in the hardcopy version of the journal (panels 'b' of Figures 1 through 13 
and 37 through 120, and Figures 13 through 36). The full data set for all metallicities is 
included in the electronic version of the journal and can be accessed on our website at 
http : //www . stsci . edu/science/starburst99/| . 



3. Spectral energy distributions 

Spectral energy distributions for all five metallicities, the three IMF parameterizations, 
and the two star-formation laws were calculated. The wavelength coverage and spectral 
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resolution are identical to those of the model atmosphere set Lejeune et al. (1997). The 
minimum and maximum wavelengths are 91 A and 160 jum, respectively. The spectral 
resolution is wavelength dependent. It is typically about 10 to 20 A in the UV to optical 
range. The energy distributions were calculated at intervals of 1, 10, and 100 Myr between 
1 and 20 Myr, 20 and 100 Myr, and 100 and 1000 Myr, respectively. 

The standard IMF case (a = 2.35; M up = 100 M Q ) is in Figure 1. Panels a-e show 
the energy distributions for the five considered metallicities. Not all time steps are plotted 
in the figures but they are included in the ascii table file which can be downloaded. The 
wavelength range has been restricted to 100 - 10,000 A in the figures, but the full range is 
covered in the tables. 

The most dramatic changes occur in the far-UV below 912 A, and even more below 
228 A in the He ++ continuum. O stars dominate during the first few Myr, followed by a 
brief period at ~4 Myr when WR stars contribute with their strong far-UV flux. Eventually 
the burst population fades away. Metallicity enters both via the continuum and line opacity 
(compare the 2600 A region at low and high metallicity) and via stellar evolution. The 
latter effect can most readily be seen in the strength of the He ++ continuum. It is strongest 
at high metallicity when WR stars preferentially form, and it is weakest at low metallicity 
(see also the discussion in § 8). 

The continuous star formation case for the same IMF is in Figure 2a-e. Most of the 
prior discussion holds for this figure as well. Once the stellar types responsible for the 
photon contribution to a wavelength interval have reached equilibrium between birth and 
death, the particular wavelength region becomes time independent. 

We move on to different IMFs. The instantaneous and continuous cases, with a = 3.3, 
M up = 100 M are in Figures 3a-e and 4a-e, respectively. The truncated IMF, a = 2.35, 
M up = 30 M is covered in Figures 5a-e and 6a-e. The spectral energy distributions 
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produce a softer radiation field as compared with the standard IMF, in particular if stars 
with M up > 30 M Q are suppressed. 

The radiative properties shown later in the paper were obtained from these energy 
distributions. For instance, colors were derived by convolving the spectra by the filter 
profiles. Users who wish to compute colors in different filter systems can do so with the 
spectra in Figures 1 through 6. Another application would be to use the spectral energy 
distributions as input to a photoionization code. In this case, Figures 1 through 6 should 
not be used because the nebular continuum was added to the stellar flux. We computed an 
additional set of stellar spectra with exactly the same parameters as before, but without 
including the nebular continuum. This set may also be useful to compute magnitudes and 
colors for comparison with objects which have no nebular continuum contribution. An 
example could be a background subtracted star cluster, where the background subtraction 
removed the nebular emission as well. There may also be cases where the nebular emission 
of a galaxy or H II region is more extended than the ionizing cluster, and part of the 
nebular flux is outside the instrument aperture. Finally, the H II region may not be 
ionization-bounded. Such cases may be easier to model using energy distributions without 
nebular continuum, and applying an empirical correction, e.g., from the observed Ha flux 
in the aperture. The models without nebular continuum are in Figure 7a-e (a = 2.35, 
M up = 100 M , instantaneous star formation), Figure 8a-e (a = 2.35, M up = 100 M , 
continuous star formation), Figure 9a-e (a = 3.3, M up = 100 M Q , instantaneous star 
formation), Figure lOa-e (a = 3.3, M up = 100 M , continuous star formation), Figure lla-e 
(a = 2.35, M up = 30 M , instantaneous star formation), and Figure 12a-e (a = 2.35, 
M up = 30 M , continuous star formation). The spectra in Figures 7 through 12 and the 
curves in Figures 45, 46, 75, and 76 are the only models in this paper which do not include 
the nebular continuum. 
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4. Ultraviolet line profiles 

Two spectral regions with useful diagnostic lines were computed at higher resolution. 
The wavelength range from 1200 A to 1850 A has the strong resonance lines of C IV A1550, 
Si IV A1400, and N V A1240. At shorter wavelengths, O VI A1035 and Ly/3 are of interest. 
The spectra described in this section were computed with a normalized library of stellar 
UV spectra and absolute fluxes derived from the energy distributions shown in the previous 
section. We did not use the energy distributions directly since they are not featureless. 
Rather we derived a featureless continuum by fitting a spline through line-free sections of 
the model atmospheres. 

4.1. The 1200 A to 1850 A region 

The 1200 A to 1850 A region is readily accessible to many satellites, like IUE or HST. 
We synthesized spectra from an IUE high-dispersion library of O- and WR-spectra and a 
low-dispersion for B stars. The spectral library and the method are discussed in Robert, 
Leitherer, & Heckman (1993) and Leitherer et al. (1995). The resolution is 0.75 A, adequate 
for typical spectral data of extragalactic objects. This resolution is only reached during 
phases dominated by O- and WR stars, as is the case for young (t < 7 Myr) bursts and for 
the continuous case at any time. Otherwise the effective spectral resolution approaches the 
resolution of IUE low-dispersion spectra (6 A) . The upgrade of the B-star library from low- 
to high- dispersion has been completed and will be implemented in the future (de Mello, 
Leitherer, & Heckman, in preparation; Robert, Leitherer, & Heckman, in preparation). 

Only solar metallicity models are shown since the library stars have solar or somewhat 
sub-solar metallicity. A spectral library of massive stars in the Magellanic Clouds is being 
developed and will be available in the future (Robert 1999; Robert, Leitherer, & Heckman 
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1999, in preparation). 

Figure 13 shows rectified spectra between 1 and 20 Myr for an instantaneous burst. 
A standard IMF was used. The evolution of the strongest stellar-wind features of C IV 
A1550, Si IV A1400, N V A1240, N IV A1720, and He II A1640 can be readily seen. These 
features disappear or become photospheric lines after about 7 Myr when the transition from 
an O-star to a B-star dominated population occurs. A purely photospheric line is, e.g., S V 
A1502. The library stars have strong interstellar lines as well. One of the strongest lines is 
C II A1335. 

The continuous case for the same IMF is shown in Figure 14. The same basic effects 
are there as well but the time-dependence is much weaker, and a quasi-equilibrium state is 
reached after about 5 Myr. Figures 15 and 16 show spectra with the same parameters as in 
Figures 13 and 14, but in absolute luminosity units. 

As we did for the spectral energy distributions in the previous section, we re-computed 
the UV spectra for the two other IMFs, a steeper IMF with a = 3.3 and a truncated IMF 
with M up = 30 M Q . These spectra are in Figures 17 through 24. The general trend is a 
weakening of all stellar-wind lines, most notably C IV A1550 and Si IV A1400, as these lines 
are due to massive O stars. These lines are very sensitive IMF tracers. Their detection 
immediately indicates an O-star population and therefore star-formation over the past 
~10 Myr. 

4.2. The 1015 A to 1060 A region 

The wavelength region shortward of Lya in local starburst galaxies is accessible to, 
e.g., the FUSE mission (Sahnow et al. 1996) and can be observed in star-forming galaxies 
at high redshift from the ground. One of the strongest stellar-wind lines is O VI A 1035 
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(Walborn & Bohlin 1996). Synthetic spectra for the region between 1015 A and 1060 A 
were presented by Gonzalez Delgado et al. (1997) and were regenerated for the Starburst99 
package. They are based on an empirical Copernicus library of O and B stars and have a 
resolution of 0.2 A. A few library stars were observed with HUT at a spectral resolution 
of about 3 A. Only solar metallicity models are considered, as we have no library stars 
with non-solar metallicity. Available time steps are 1 to 20 Myr at At = 1 Myr and 20 to 
100 Myr at At = 10 Myr. Not every time step is plotted in the figures. 

We discuss the most notable features using the standard IMF, instantaneous burst, 
rectified model (Figure 25). The strongest lines are O VI A1035, Ly/3, C II A1036, and 
the Lyman and Werner bands of H2, like the one at 1050 A. In the first few Myr of the 
starburst, O VI A1035 has a strong P Cygni profile, indicating winds from massive O stars. 
O VI A1035 becomes weaker with age and gradually blends with Ly/3 at 1026 A. Ly/3 is 
mostly stellar, with some interstellar contribution. The line increases in strength with 
increasing B-star fraction. C II A 1036 follows Ly/5 in its temporal behavior. It is a strong 
B-star line with some interstellar contribution and can be used to estimate the starburst 
age (Gonzalez Delgado et al. 1997). The continuous case with the same IMF parameters 
behaves accordingly (Figure 26). The counterparts of Figures 25 and 26 in luminosity units 
are in Figures 27 and 28. 

The remaining figures for O VI A1035 are organized the same way as the figures for the 
1200 A to 1850 A region. They are intended to highlight IMF variations (Figures 29 through 
36). An IMF biased against massive O stars suppresses O VI A1035 and strengthens C II 
A1036. Even at the earliest ages the absorption trough around 1025 A is mostly due to 
Ly/3, and not O VI A1035. 
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5. Massive-star inventory 

In this section we give number predictions for massive stellar types that can be easily 
"counted" even in distant galaxies: stars, WR stars, and supernovae (SNe). 

O-star numbers are in Figures 37a-e and 38a-e for the instantaneous and continuous 
case, respectively Panels a, b, c, d, and e are for metallicities 0.040, 0.020, 0.008, 0.004, 
and 0.001, respectively. All remaining figures in this paper have the same structure. The 
O-star numbers include stars with spectral types 03 to 09.5 of all luminosity classes. The 
adopted spectral-type versus effective temperature (T c g) and luminosity (L) relation is from 
Schmidt-Kaler (1982). The chosen star-formation histories produce up to about 10 3 to 10 4 
O stars at any time, depending on the age of the population. A flatter IMF increases the 
O-star number. 

WR stars are the evolved descendants of massive O stars. We define them as stars with 
logT c fj > 4.4 and surface hydrogen abundance less than 0.4 by mass. The mass limits for 
the formation of WR stars are taken from Maeder & Meynet (1994). WR stars have very 
extended atmospheres producing strong emission lines which can be detected in distant 
starburst populations (Conti 1991). The quantity of interest for comparison with models 
is the WR/O ratio, which is plotted in Figures 39a-e and 40a-e for the instantaneous 
and continuous case, respectively. This ratio is very metallicity dependent: metal-rich 
starbursts are predicted to have more WR stars at any time and to have longer phases 
when WR stars are present. This behavior is generally in agreement with observations 
(Meynet 1995). The panels with sub-solar metallicity have no graph for the IMF with 
M up = 30 M . This results from the absence of WR stars at low metallicity. The detailed 
shape of the graphs in Figures 39a-e and 40a-e is very model dependent and future revisions 
of the evolutionary models may change the WR/O ratio. Dedicated evolutionary synthesis 
models with particular attention to the WR phase were published by Meynet (1995) and 
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Schaerer & Vacca (1998). The binary channel to form WR stars in stellar populations was 
investigated by Cervino & Mas-Hesse (1997), Vanbeveren et al. (1997), Schaerer & Vacca 
(1998), and Dionne (1999). 

WR stars come in two main subclasses: nitrogen-rich WN stars and carbon-rich WC 
stars. We adopt the classification schemes of Conti, Leep, & Perry (1983) for WN stars and 
of Smith & Hummer (1988) and Smith & Maeder (1991) for WC stars. The WC/WN ratio 
measures the relative lifetimes spent in the two phases (Figures 41a-e and 42a-e). However, 
the calculations are affected by uncertainties in the interpolation of the relatively sparse 
tracks, as discussed in Schaerer & Vacca (1998). In this respect, the present models are 
identical to those of Schaerer & Vacca, who predict fewer WC stars than Meynet (1995) 
despite the use of the same tracks. Most of the comments on the WR/O ratio apply to 
these figures as well. The WC/WN ratio varies very little with the IMF exponent: the 
a = 2.35 and 3.3 cases are almost indiscernible in the figures. This indicates that both 
stellar types have progenitors of very similar mass. As a result, applying different weights 
to different mass intervals does not affect the ratio. The same argument applies to several 
other quantities discussed further below, like, e.g., the CO index (Figure 95). 

All massive stars with initial masses above 8 M are assumed to explode as SNe. This 
leaves open the question if a critical mass exists above which stars directly form a black 
hole rather than producing a core-collapse SN (Maeder 1992). If so, the predictions of 
Figures 43a-e and 44a-e would hardly be changed because the SN rate is strongly weighted 
towards the lowest progenitor masses. 

The user should be aware of a computational issue which can be seen in Figures 43a-e 
and 44a-e. The graphs show discontinuities at the factor-of-2 level where there should 
be none (e.g., around 10 Myr in Figure 43). A simple argument suggests that the SN 
rate should be a smooth and almost constant function of time. The rate scales with the 
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product of the IMF and the mass-age relation. For a Salpeter IMF and a mass-age relation 
t(m) oc wT 1 with 7 rs 1.2 (Schaerer et al. 1993b) the SN rate of an instantaneous burst 
becomes essentially time-independent (Shull & Saken 1995). In other words, there are more 
and more SN events for lower masses since there are more progenitors available, but at 
the same time it takes longer for a star to turn into a SN. Obviously the discontinuities 
in Figure 43 are not physical. They result from a peculiarity of the isochrone synthesis 
interpolation technique. The code attempts to optimize the mass interval size for the 
interpolation. This works extremely well, except when only stars in an infinitesimal mass 
interval are relevant, as is the case for the SN rate. The same applies to all quantities which 
are directly dependent on the SN rate, like the energy release from SNe. We decided not to 
artificially smooth the curves since this would be a subjective process and we did not want 
to treat SN-related quantities different from the other predictions. We remind the user to 
apply common sense when interpreting the figures. The correct supernova rate between 5 
and 35 Myr for an instantaneous burst in Figure 43 is about 10~ 3 yr _1 (LH95). 

6. Luminosities 

We give the total and monochromatic luminosities at selected wavelengths. We have 
chosen the V, B, and K passbands, as well as a UV wavelength at 1500 A since these are 
the most interesting cases for comparison with observations. Luminosities for other bands 
can be obtained from the spectral energy distributions in Figures 1 through 6. 

The bolometric luminosity (Mb i) is defined such that the Sun has Mb i = 4.75. 
The results in Figures 45a-e and 46a-e were obtained by integrating the spectral energy 
distributions without the nebular continuum. The total luminosity of the combined stellar 
and nebular continuum is smaller than the pure stellar continuum since only a fraction of 
the absorbed stellar ionizing flux is re-emitted in the nebular continuum. The curves are 



-19- 



smooth without strong discontinuities and reflect the behavior of the most massive stars 
which provide most of the radiative energy output. 

The absolute magnitudes My, Mb, and Mk are in Figures 47a-e (instantaneous) 
through 52a-e (continuous). The absolute magnitude My was calculated from the 
bolometric luminosity and the bolometric correction, which is 0.00 for a solar metallicity 
star with T c g = 7000 K and log g = 1.0. The Sun has a bolometric correction of —0.19 
in this system. The absolute magnitudes Mb and Mk follow from My and the (B — V) 
and (V — K) colors. All three luminosities have conspicuous variations around 10 Myr 
when red supergiants (RSG) appear. The effect is strongest in the K band, which is 
closest to the energy peak of RSGs. The RSG feature is very metallicity-dependent. The 
RSG contribution in Figure 51a-e is strongest at Z — 0.040 and weakest at Z — 0.001. 
This apparent metallicity dependence results from the failure of the evolutionary models 
to predict correct RSG properties at Z < 0.008. This effect was pointed out before by 
Mayya (1997) and Origlia et al. (1998). Low-metallicity RSG models have too high surface 
temperatures and too short RSG lifetimes. Uncertainties in the mixing processes at low 
metallicity are a possible explanation (Langer & Maeder 1995). Currently there are no 
self-consistent stellar evolution models available which correctly predict the variations of 
blue-to-red supergiants with metallicity. Therefore our (and most other) synthesis models 
are incorrect during phases when RSGs are important. An empirically adjusted set of 
synthesis models was prepared by Origlia et al. (1998) but these adjustments were not made 
for the model set in this paper. 

Luminosities at 1500 A (Figures 53a-e and 54a-e) were computed by averaging over 
the wavelength interval 1490 A - 1510 A. This wavelength becomes observable from the 
ground at redshifts larger than ~2, and the luminosity at 1500 A is a useful indicator of the 
star-formation rate, or, in connection with the Ha luminosity of dust obscuration (Pettini 
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et al. 1997). The 1500 A luminosity is a very robust prediction, with few uncertainties since 
the continuum at this wavelength comes from well-understood late-O/early-B stars. 

7. Colors 

Optical and near-infrared (IR) colors were calculated by convolving the spectral energy 
distributions with the filter profiles. The filters are in the Johnson (1966) system and are 
the same as in LH95. The zero point is defined by a star with Z = 0.020, T c g = 9400 K, 
and log g = 3.95, which has zero colors in all passbands. Colors in other photometric 
systems can be obtained by convolving the spectral energy distributions with the desired 
filter profiles. 

These colors are available: (U — B) (Figures 55a-e and 56a-e), (B — V) (Figures 57a-e 
and 58a-e), {V — R) (Figures 59a-e and 60a-e), (V — I) (Figures 61a-e and 62a-e), (V — J) 
(Figures 63a-e and 64a-e), (V — H) (Figures 65a-e and 66a-e), (V — K) (Figures 67a-e 
and 68a-e), and (V — L) (Figures 69a-e and 70a-e). As a reminder, the continuous nebular 
emission is included in the colors, but not the line emission. The R filter is the most likely 
passband to suffer from nebular line contamination since Ha is included. A first check of 
the expected degree of contamination can be made by comparing the Ha equivalent widths 
of Figures 83a-e and 84a-e with the width of the R filter (about 2000 A). Line strengths 
of other strong lines in H II regions can be estimated from the photoionization models of 
Stasihska & Leitherer (1996). 

The RSG issue raised before applies to some of the color plots for low metallicities as 
well. The strong metallicity dependence of the RSG feature around 10 7 yr is related to the 
effects discussed in § 6. 

Rather than colors, we give continuum slopes in the UV. We define the slope /3 as 



-21- 

the spectral index of the spectral energy distribution: F\ oc X 13 . Two indices are shown, 
one for the average slope between 1300 A and 1800 A (Figures 71a-e and 72a-e), and 
one for the 2200 A to 2800 A region (Figures 73a-e and 74a-e). The slopes were simply 
derived by fitting a first-order polynomial to the spectra through the wavelength intervals 
1280 A - 1320 A and 1780 A -1820 A for /3 1550 , and through 2180 A - 2220 A and 
2780 A - 2820 A for /^soo • This should serve as an approximate guide for the variation of 
the slope with time but becomes increasingly meaningless if the actual spectrum deviates 
from a power law. Figures 1 through 6 suggest that spectral energy distributions with ages 
less than ~200 Myr are indeed well approximated by a power law in the UV but that this 
assumption is no longer correct at older ages. The UV slopes are quite independent of 
evolution and IMF effects for the first 30 Myr. This property makes them very useful for 
deriving UV extinctions in galaxy spectra (Calzetti, Kinney, & Storchi-Bergmann 1994). 
Therefore the spectral slope of the restframe ultraviolet spectra of star-forming galaxies at 
high redshift can be used to estimate the effects of dust obscuration in the early universe 
(Calzetti & Heckman 1999, and references therein). 

8. Far-ultraviolet properties 

The predictions in this section rely on our capabilities to model the stellar far-UV 
continuum below 912 A since this spectral region is generally not accessible to direct 
observations. A comparison between different models for hot stars in this spectral region 
has been made by Schaerer & de Koter (1997). The internal consistency, as judged from 
differences between the models, is within 0.1 dex between 912 A and 504 A, which is 
relevant for the flux ionizing neutral hydrogen (JV(H°)). Uncertainties become larger 
towards shorter wavelengths, in particular below 228 A, where the emergent flux becomes 
strongly dependent on stellar-wind properties. This of course leaves open the question of 
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the external, absolute uncertainties. Model atmospheres generally do a good job above 
228 A when combined with photoionization models and compared to H II region spectra 
(Garcia- Vargas 1996) although there are some properties which are sensitive to the flux 
distribution in the neutral He continuum (Stasihska & Schaerer 1997). This makes errors 
in the integrated photon fluxes by more than 0.3 dex unlikely. The region below 228 A has 
not yet been tested in such detail and the uncertainties are potentially large. 

Far-UV models for stars colder than about 30,000 K (types B and later) are lagging 
behind. Cassinelli et al. (1995) found discrepancies by two orders of magnitude between 
the observed extreme-UV flux of the B2 star e CMa and model predictions. Even slight 
errors in the adopted wind parameters can produce huge far-UV flux variations at these 
relatively low temperatures (Schaerer & de Koter 1997). 

We define the Lyman break as the ratio of the average flux in the wavelength interval 
1080 A - 1120 A over that between 870 A and 900 A. These wavelength intervals are 
relatively free of line-blanketing so that we measure mostly temperature rather than line 
opacity. The Lyman break is about a factor of 2 to 3 in O-star dominated phases and 
drops thereafter (Figures 75a-e and 76a-e). The warnings about model uncertainties for the 
Lyman continuum of cooler stars apply to the instantaneous case. Figure 75a-e should not 
be overinterpreted after about 50 Myr since B stars like e CMa (see previous paragraph) 
could dominate. The Lyman break for a continuous population is always dominated by hot 
stars so that Figure 76a-e can be trusted over the entire range plotted. 

The number of photons capable of ionizing neutral hydrogen (JV(H°)), neutral helium 
(JV(He°)), and ionized helium (iV(He + )) were calculated by integration of the spectra below 
912 A, 504 A, and 228 A, respectively. They are shown in Figures 77a-e and 78a-e (JV(H°)), 
Figures 79a-e and 80a-e (7V(He°)), and Figures 81a-e and 82a-e (iV(He + )). As stated before, 
the numbers become increasingly uncertain with shorter wavelength. iV(He + ) is almost 
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entirely produced by WR stars, and uncertainties in the wind properties enter. 

We also give equivalent widths of several popular hydrogen recombination lines. The 
continuum is taken from the model atmospheres and does not include underlying stellar 
absorption. The nebular continuum is of course taken into account. The individual plots 
are for Ha (Figures 83a-e and 84a-e), H/3 (Figures 85a-e and 86a-e), Pa/5 (Figures 87a-e 
and 88a-e), and Br7 (Figures 89a-e and 90a-e). The transformation relations from iV(H°) 
to the line luminosities are in LH95. 



9. Other diagnostic lines 

In this section we present a few more diagnostics that can be useful to isolate a stellar 
population in a galaxy spectrum. They are all related to stars off the main-sequence: WR 
stars, RSGs and SNe. 

In Figures 81 and 82 we predicted iV(He + ), which can immediately be converted into 
the emission-line flux of nebular He II A4686. Generally, a hot-star population capable of 
producing nebular He II will also show broad, stellar He II A4686 (Schaerer & Vacca 1998). 
Our model predictions for the stellar line are in Figures 91a-e and 92a-e. The feature is 
commonly referred to as the "WR bump" . It includes only He II and none of the other 
nearby spectral features like C III, N III, [Ar IV], and [Fe III]. Note that our code predicts 
other WR lines from the list of Schaerer & Vacca but they are not given here. 

[Fe II] A1.26 is useful to count supernovae in starbursts (Figures 93a-e and 94a-e). The 
supernova shock wave destroys interstellar dust grains, thereby releasing iron atoms and 
ions which had condensed on the dust grains. We adopted the scaling relation of Calzetti 
(1997) to convert the supernova rates into [Fe II] line luminosities. 

The remaining figures in this section are related to RSG properties. We repeat again 
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that the predicted properties of post-main-sequence stars are not reliable and that all 
phases dependent on RSG properties at sub-solar metallicity are suspect. We begin with 
the spectroscopic CO index at 2.2 /iin (Figures 95a-e and 96a-e). We follow the definition 
of Doyon, Joseph, & Wright (1994) who expressed the CO strength as a function of 
temperature for dwarfs, giants, and supergiants. The index is set to for T c g > 6000 K. 
We also computed the strength of the calcium triplet at AA8498 A, 8542 A, 8662 A 
(Figures 97a-e and 98a-e). The equivalent width of the sum of all three lines is related to 
gravity (log g) and metallicity by the relation W(CaT) = 10.21 — 0.95 log g + 2.18 log Z/Z & 
(Diaz, Terlevich, & Terlevich 1989; Garcia- Vargas, Molla, & Bressan 1998). The predicted 
values neither take into account nebular emission nor stellar absorption of higher Paschen 
lines which fall in this wavelength region. Observations must be corrected for these lines, if 
present, before a comparison is made. 

Evolutionary synthesis models for the first and second overtones of CO at 2.29 /xm and 
1.62 /zm and for Si I at 1.59 /xm were presented by Origlia et al. (1998). The models in 
Figures 99a-e and lOOa-e (CO A1.62 /zm), Figures lOla-e and 102a-e (CO A2.2 /im), and 
in Figures 103a-e and 104a-e (Si I A1.59 /zm) are expanded versions of the unmodified case 
discussed by Origlia et al. In that paper, the effect of modifying the evolutionary tracks 
was studied and several sets of IR-line models were published. To be consistent with the 
other quantities shown here, we opted for including only the standard models. The models 
for Si I A1.59 /im were not covered by Origlia et al. and are shown here for the first time. 
The theoretical library of Origlia et al. (1993) which was used for Si I is less reliable than 
that for CO so that care is required when using the Si I models. 
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10. Mass and energy return 

The previous sections cover stellar numbers and radiative properties. Here we turn to 
non-radiative properties of the stellar population. The input physics is discussed in greater 
detail in Leitherer et al. (1992). The figures in this section show the mass and energy input 
by stars and SNe. Only core-collapse SNe are considered. They are assumed to release 
10 51 erg per event in the form of kinetic energy, independent of metallicity. We do not 
address the efficiency of thermalization, which would require hydrodynamical modeling. 
The simulations of Thornton et al. (1998) suggest that about 10% of the available kinetic 
energy can actually be used to pressurize the interstellar gas. The remaining 90% are 
radiated away. 

The rate of mass return of stellar-wind and SN material is plotted in Figures 105a-e 
and 106a-e. The adopted mass-loss rates are not those of the evolutionary models but 
those favored by Leitherer et al. (1992) and LH95. We prefer this approach over simply 
using the evolutionary mass-loss rates. In our opinion, the evolutionary mass-loss rates 
are too high (see Section 2) and should be considered only as an adjustable parameter 
in evolution models but arc not directly related to the observed mass-loss rate. Slight 
differences between Figures 105 and 106 and the corresponding figures in LH95 are not due 
to a different mass-loss parameterization but because of different stellar parameters (T e g, L) 
in the revised tracks. The individual contributions from stellar winds and SNe are broken 
down in Figures 107a-e and 108a-e for the standard IMF case. Stellar winds are generally 
more important for young bursts whereas SNe take over at later times. 

The total mass return from winds and SNe is in Figures 109a-e and llOa-e. The 
quantity plotted is / Molt, i.e. the integral of the curves in Figures 107 and 108 over time. 
This quantity is useful to evaluate the exhaustion of the gas supply in a galaxy or the 
degree of chemical pollution by wind and supernova material. 
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Figures llla-e and 112a-e give the mechanical luminosity L mec h released by winds and 
supernovae. The curves are similar to those for the mass return in Figures 105 and 106. The 
relative contributions to L mec h from stellar winds and SNe are in Figures 113a-e and 114a-e. 
A further break-down into individual stellar-wind components is given in Leitherer et a!. 
(1992). Generally, most of the wind power comes from WR stars, with some contribution 
from O stars. All other stellar phases are negligible since wind velocities of cool stars are 
lower by two orders of magnitude. The energy return (/ L mec h^) from winds and SNe is in 
Figures 115a-e and 116a-e. 

It is instructive to perform a differential comparison between the radiative and the 
non-radiative energy output of young stellar populations. This is done in Figures 117a-e 
and 118a-e for the ratio of the ionizing (<912 A) over the bolometric luminosity, and 
in Figures 119a-e and 120a-e for the ratio of the mechanical (L mec h) over the bolometric 
luminosity. The non-radiative energy input into the interstellar medium becomes significant 
in comparison with ionizing radiation once the strong WR winds are turned on. For a 
single population, non-thermal and ionizing energy input become equally important around 
~10 Myr. 

11. Conclusions 

We have computed a large grid of predictions for observable properties of galaxies 
with active star formation. The distribution of the models is purely web-based. We believe 
the community will find this method more useful than a hardcopy publication. All the 
figures discussed in this paper are at pittp : //www . stsci . edu/science/starburst99/| . 
This webpage provides links to other spectrophotometric databases as well. 

It is worthwhile to recall the most important short-comings and uncertainties of the 
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models: 

Chemical evolution is not treated self-consistently. Each stellar generation has the same 

chemical composition. This becomes a concern for models which are evolved over times 

during which significant changes in the metallicity of the ISM occur. In this case spectra 

will have a wavelength dependent metallicity. Generally, light at shorter wavelengths is 

produced by more massive and younger stars, which are chemically more evolved than older 

stars. 

Binary evolution has been neglected. Although about 50% of all stars form in binaries, it is 

not clear if a significant fraction of these binaries has an evolutionary history that differs 

substantially from single-star evolution. There is disagreement in the literature on this 

point. Under these circumstances we took the approach of implementing the simpler of the 

two alternatives, as nature itself often prefers to do. 

Mass loss and mixing processes in stellar evolution are still poorly understood. Stellar 

phases, like WR stars or RSGs, are particularly affected by such uncertainties. The main 

culprit is the lack of a self-consistent theory which makes the introduction of adjustable 

parameters necessary. These parameters then are sometimes extrapolated into a regime for 

which they were not calibrated, such as metallicity. 

The situation is particularly disturbing for red supergiants whose properties are badly 

reproduced at low metallicity. There is no easy work-around for the user of evolutionary 

models, except for an empirical adjustment of the tracks. We have discussed such an 

adjusted model set. Clearly, a strong effort on the stellar evolution modeling side is called 

upon for improvements. 

Our models put most of the emphasis on early evolution phases. Later phases, like AGB 

stars or white dwarfs are covered only crudely or not at all. While improvements to our 

code can be made, our prime goal and expertise is related to massive, hot stars, and we 

decided to optimize this stellar species first. 
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Despite the warnings, the model set should turn out to be useful for the interpretation 
of observations of star- forming galaxies. For maximum benefit, the user is encouraged to 
compare our model predictions with those of other groups, such as those mentioned earlier 
in the paper. 

To provide maximum flexibility, we offer the user to run tailored models at our website. 
Instructions on how to run the code and how to access the results are given at the website. 
The Fortran code is distributed freely and can be retrieved from the website as well. 

Harry Payne helped us trouble-shoot numerous bugs and pitfalls we encountered 
during the construction of our webpage. D. Foo Kune acknowledges financial support from 
the STScI Summer Student Program. Salary support for D. Devost and D. Schaerer came 
from the STScI Director's Discretionary Research Fund. We appreciate advice on electronic 
publishing and data maintenance from Bob Hanisch. Partial support for this work was 
provided by NASA through grant number NAG5-6903, from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc., under NASA contract NAS5-26555. 



-29- 

REFERENCES 

Calzetti, D. 1997, AJ, 113, 162 

Calzetti, D., & Heckman, T. M. 1999, ApJ, in press 

Calzetti, D., Kinney, A. L., & Storchi-Bergmann, T. 1994, ApJ, 429, 582 

Carico, D. P., Sanders, D. B., Soifer, B. T., Elias, J. H., Mathews, K., & Neugebauer, G. 
1988, AJ, 95, 356 

Cassinelli, J. P., Cohen, D. H., MacFarlane, J. J., Drew, J. E., L A. E., Hoare, M. C, 

Vallerga, J. V., Welsh, B. Y., Vedder, P. W., Hubeny, I., & Lanz, T. 1995, ApJ, 438, 
932 

Cervifio, M., & Mas-Hesse, J. M. 1994, A&A, 284, 749 

Cervino, M., Mas- Hesse, J. M., & Kunth, D. 1997 in Starburst Activity in Galaxies, ed. J. 
Franco, R. Terlevich, & A. Serrano, Rev. Mex. Astron. Astrofis. Conf. Ser., 6, 188 

Charbonnel, D., Meynet, G., Maeder, A., Schaller, G., & Schaerer, D. 1993, A&AS, 101, 
415 

Chariot, S., & Bruzual G. A. 1991, ApJ, 367, 126 

Conti, P, S. 1991, ApJ, 377, 115 

Conti, P. S., Leep, E. M., & Perry, D. N. 1983, ApJ, 268, 228 

Diaz, A. L, Terlevich, E. & Terlevich, R. 1989, MNRAS, 239, 325 

Dionne, D. 1999, in IAU Symp. 193, Wolf-Rayet Phenomena in Massive Stars and Starburst 
Galaxies, ed K. A. van der Hucht, G. Koenigsberger, & P. R. J. Eenens (San 
Francisco: ASP), in press 

Doyon, R., Joseph, R. D., & Wright, G. S. 1994, ApJ, 421, 101 

Ferland, G. J. 1980, PASP, 92, 596 



-30- 

Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 

Garcia- Vargas, M. L. 1996, in From Stars to Galaxies: The Impact of Stellar Physics on 
Galaxy Evolution, ed. C. Leitherer, U. Fritze-von Alvensleben, & J. Huchra (San 
Francisco: ASP), 244 

Garcia- Vargas, M. L., Molla, M., & Bressan, A. 1998, A&AS, 130, 513 

Girardi, L., & Bertelli, G. 1998, MNRAS, in press 

Gonzalez Delgado, R. M., Leitherer, C, & Heckman, T. M. 1997, ApJ, 489, 601 

Johnson, H. L. 1966, ARAA, 4, 193 

Kurucz, R. L. 1992, in IAU Symp. 149, The Stellar Populations of Galaxies, ed. B. Barbuy 
& A. Renzini (Dordrecht: Kluwer), 225 

Langer, N., & Maeder, A. 1995, A&A, 295, 685 

Leitherer, C. 1998, in 38th Herstmonceux Conference, The Stellar Initial Mass Function, 
ed. G. Gilmore, I. Parry, & S. Ryan (San Francisco: ASP), 61 

Leitherer, C, et al. 1996a, PASP, 108, 996 

Leitherer, C, Chapman, J., & Koribalski, B. 1995, ApJ, 481, 898 

Leitherer, C, Fritze-v. Alvensleben, U., & Huchra, J. 1996b, From Stars to Galaxies -- The 
Impact of Stellar Physics on Galaxy Evolution (San Francisco: ASP) 

Leitherer, C, & Heckman, T. M. 1995, ApJS, 96, 9 (LH95) 

Leitherer, C, Robert, C, & Drissen, L. 1992, ApJ, 401, 596 

Leitherer, C, Robert, C, & Heckman, T. M. 1995, ApJS, 99, 173 

Lejeune, T., Buser, R., & Cuisinier, F. 1997, A&AS, 125, 229 

Maeder, A. 1990, A&AS, 84, 139 

-. 1992, A&A, 264, 105 



-31- 

— . 1995, in Astrophysical Applications of Stellar Pulsation, ed. R. S. Stobie & P. A. 
Whitelock (San Francisco: ASP), 1 

Maeder, A., & Meynet, G. 1994, A&A, 287, 803 

Mayya, Y. D. 1997, ApJ, 482, 149 

Meynet, G. 1995, A&A, 298, 767 

Meynet, G., Maeder, A., Schaller, G., Schaerer, D., & Charbonnel, C. 1994, A&AS, 103, 97 

Moller, C, Fritze-v. Alvensleben, U., & Fricke, K. J. 1997, A&A, 317, 676 

Moorwood, A. F. M. 1996, Space Sci. Rev., 77, 303 

Origlia, L., Goldader, J. D., Leitherer, C, Schaerer, D., & Oliva, E. 1998, ApJ, in press 

Origlia, L., Moorwood, A. F. M., & Oliva, E. 1993, A&A, 280, 536 

Pettini, M., Steidel, C C, Dickinson, M. E., Kellog, M., Giavalisco, M., Adelberger, K. L. 
1997, in The Ultraviolet Universe at Low and High Redshift, ed. W. H. Waller, M. 
N. Fanelli, J. E. Hollos, & A. C. Danks (Woodbury: AIP), 279 

Robert, C. 1999, in IAU Symp. 193, Wolf-Rayet Phenomena in Massive Stars and Starburst 
Galaxies, ed K. A. van der Hucht, G. Koenigsberger, & P. R. J. Eenens (San 
Francisco: ASP), in press 

Robert, C, Leitherer, C, & Heckman, T. M. 1993, ApJ, 418, 749 

Sahnow, D. J., Friedman, S. D., Oegerle, W. R., Moos, H. W., Green, J. C, & Siegmund, 
O. H. 1996, Proc. SPIE, 2807, 2 

Salpeter, E. E. 1955, ApJ, 121, 161 

Scalo, J. 1998, in 38th Herstmonceux Conference, The Stellar Initial Mass Function, ed. G. 
Gilmore, I. Parry, & S. Ryan (San Francisco: ASP), 201 

Schaerer, D., Charbonnel, C, Meynet, G., Maeder, A., & Schaller, G. 1993a, A&AS, 102, 
339 



-32- 

Schaerer, D., & de Koter, A. 1997, A&A, 322, 598 

Schaerer, D., Meynet, G., Maeder, A., & Schaller, G. 1993b, A&AS, 98, 523 

Schaerer, D., & Vacca, W. D. 1998, ApJ, 497, 618 

Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 96, 269 

Schmidt-Kaler, T. 1982, in Landolt-Bornstein, New Series, Group VI, Vol. 2b, ed. K. 
Schaifers & H. H. Voigt (Berlin: Springer), 1 

Schmutz, W., Leitherer, C, & Gruenwald, R. 1992, PASP, 104, 1164 

Shull, J. M., & Saken, J. M. 1995, ApJ, 444, 663 

Smith, L. F., & Hummer, D. G. 1988, MNRAS, 230, 511 

Smith, L. F., & Maeder, A., 1991, A&A, 241, 77 

Stasihska, G., & Leitherer, C. 1996, ApJS, 107,427 

Stasihska, G., & Schaerer, D. 1997, A&A, 322, 615 

Thornton, K., Gaudlitz, M., Janka, H.-Th., & Steinmetz, M. 1998, ApJ, 500, 95 

Vanbeveren, D., Van Bever, J., & de Donder, E. 1997, A&A, 317, 487 

Walborn, N. R., & Bohlin, R. C. 1996, PASP, 108, 477 

Weedman, D. W. 1987, in Star Formation in Galaxies, ed. C. J. Lonsdale (Washington: 

NASA), 351 

Weedman, D. W., Feldman, F. R., Balzano, V. A., Ramsey, L. W., Sramek, R. A., & Wu, 
C.-C. 1981, ApJ, 248, 105 



This manuscript was prepared with the AAS IATjtN. macros v4.0. 



